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A reversed- field  plasma  confinement  geometry  has  been  produced  by  plasma  currents 
induced  by  a rotating  relativistic  electron  beam  injected  into  85  mTorr  neutral  hydrogen.  The 
resulting  belt-pinch-like  configuration  has  a rectangular  cross-section,  with  length-to-wldth  ratio 
of  20,  and  a midplane  p ~ .5.  An  equilibrium  model,  fitted  to  measurements  of  the  excluded 
flux  and  induced  axial  (JiB,)  and  azimuthal  (Bj;)  magnetic  fields,  is  in  agreement  with  the  annular 
plasma  profile  observed  with  radial  Thomson  scattering  scans,  CO2  interferometry  and  end-on  streak 
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GENERATION  OF  A REVERSED-FIELD  PLASMA 
CONFINEMENT  CONFIGURATION  WITH  A ROTATING  RELATIVISTIC 

ELECTRON  BEAM 


I.  Introduction 

It  has  recently  been  proposed  to  produce  a reversed-field  plasma 
confinement  configuration  with  a rotating  relativistic  electron  beam.^ 
This  paper  describes  initial  studies  of  such  a configuration,  which  is 
similar  to  an  elongated  cross-section  Tokamak  or  a Belt  Pinch^  (see 
Fig.  1) , and  is  formed  by  beam-induced  azimuthal  and  axial  plasma 
currents  of  sufficient  magnitude  to  reverse  the  applied  axial  (poloidal) 
field  on  axis,  and  produce  a comparable  azimuthal  (toroidal)  field, 
respectively.  As  the  induced  currents  persist  long  after  beam  passage, 
a plasma  confined  by  such  a configuration,  once  set  up,  can  be  heated 
further  by  magnetic  compression  (for  example,  by  imploding  metallic 
liners^) , or  by  injection  of  high  energy  neutral  beams  (as  in  2 x IIB"*)  . 

The  formation  of  this  reversed-field  configuration  is  shown 
schematically  in  Figs.  2a-2d.  (For  simplicity,  only  the  aximuthal 
currents  and  axial  magnetic  fields  are  shown.)  A rotating  relativistic 
electron  beam,  produced  by  passing  an  annular  beam  through  a magnetic 
cusp,  is  injected  along  an  applied  magnetic  field,  Bg,  into  a conduct- 
ing chamber  containing  neutral  hydrogen  gas.  The  large  dB/dt  at  the 
Note;  Manuscript  submitted  November  9,  1977. 
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beam  front  Induces  an  electric  field  that  breaks  down  the  gas  and  induces 

a return  current,  i.e.  plasma  currents  in  a direction  counter  to  the 

beam  current.  The  neutral  gas  pressure  is  chosen  such  that  breakdown 

is  fast  enough  to  allow  a sufficiently  dense  plasma  for  beam  propagation, 

but  slow  enough  to  keep  the  rise  in  conductivity,  and  consequently  the 

return  current,  small.  In  other  words,  the. beam  is  charge*,  but  not 

current-,  neutralized.  Thus, large  net  current  densities  will  be 

generated  across  along  the  applied  magnetic  field, 

inducing  an  axial  (AB  ) and  an  azimuthal  (B.)  magnetic  field, 

2 0 

respectively  (Fig.  2a).  The  relative  magnitudes  of  J.  and  J 

will  depend  on  the  ratio  of  the  azimuthal  (v^)  and  axial  (v^) 
velocity  components  of  the  beam  electrons.  During  the  beam  pulse,  the 
plasma  is  heated  by  direct  beam-plasma  interactions^  as  well  as  ohmic  dis- 
sipation of  plasma  currents^  (Fig.  2b).  As  the  beam  current  decreases, 

Che  dB/dc  of  the  decaying  beam  induces  a change  in  plasma  current  that 
opposes  the  change  in  beam  current.  The  high  conductivity  of  the  heated 
plasma  demands  chat  these  changes  be  almost  equal,  in  order  to  conserve 
magnetic  flux  (Fig.  2c)  , thus  maintaining  the  current  configuration  estab- 
lished during  the  beam  pulse.  For  a sufficiently  large 
AB^  is  strong  enough  Co  produce  reversal  of  Che  applied  field  on  axis 

(i.e.  IaB  I 2 1.0  B ).  If  Che  chamber  has  conducting  end  walls,  Che 
* z * o 

field  lines  must  close  inside  them,  and  Che  reversed- field  geometry 

shown  in  Fig.  2d  is  formed.  The  magnitude  of  AB^  can  be  calculated 

from  the  net  axial  current,  I the  azimuthal  to  axial  beam 

z-net 

velocity  ratio,  v./v  , and  the  layer  radius,  r~; 
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The  comb inac Ion  of  che  applied  axial  field,  che  induced  axial  field, 
and  che  induced  azimuchal  field,  leads  co  an  elongated  crosS'SecCion 
toroidal  plasma  configuration.  Such  a field  configuration  is  expected 
to  exhibit  favorable  confinement  properties  for  high  9 plasmas.^ 

The  condition  for  the  sequence  of  events  described  above  Co  occur  is; 


”^01  "^beam  ”^02’ 


where  t. 


is  che  beam  duration,  and  t-,  and  t _ are  che  magnetic 


field  decay  times  in  Che  plasma  before  and  after  beam  passage,  respectively. 
For  diffusive  decay,  is  related  to  the  plasma  conductivity,  j^, 

and  che  characteristic  length,  x,  by 


^D1  ' ‘“<>‘’1’'“' 


It  can  be  seen  that  for  a typical  relativistic  electron  beam  of  5O  nsec 
duration,  injected  into  a system  with  a characteristic  length  of  a few 
cm,  che  target  gas  must  be  Initially  neutral  or  only  very  weakly 
ionized  to  satisfy  the  first  inequality  in  Eq.  (2),  whereas  a plasma 
with  final  electron  temperature  exceeding  1 eV  is  sufficient  to  satisfy 
che  second.  Thus,  che  beam  is  not  required  to  deposit  much  of  its 
energy  as  heat,  but  rather  to  act  as  a vehicle  to  transport  magnetic 
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energy  Into  che  system,  which  In  turn  remains  there  long  after  beam 
passage . 

The  Induced  plasma  currents  will  decay  on  a time  scale  determined 
by  the  plasma  conductivity.  Therefore,  in  the  absence  of  significant 
loss  mechanisms,  their  dissipation  will  result  in  further  plasma  heating, 
and  che  configuration  should  persist  for  Che  characteristic  magnetic 
diffusion  time  of  the  geometry.  It  can  be  seen,  Chen, that  in  this  system 
the  induced  magnetic  field  plays  a dual  role;  1)  it  provides  che  plasma 
confinement  and  2)  it  acts  as  a store  for  che  beam's  energy,  which  can 
then  be  converted  into  plasma  thermal  energy  through  resistive  dissipation. 
This  latter  effect  was  proposed  and  Chen  observed  in  earlier  rotating 
beam  experiments  by  KapeCanakos,  ^ which  achieved  field  reversed 

configurations  having  a duration  of  less  chan  I usee. 

In  Che  rotating  beam  experiment  of  Roberson,  ^ strong 

reversals  of  che  applied  field  (|aB^|  > with  lifetimes  of  2 usec 

were  observed.  It  was  also  found  that  I _ v,/v  , and  consequently 
AB^,  was  independent  of  Che  applied  magnetic  field. 

To  summarire  che  experimental  results  presented  here,  che  formation 
of  an  annular  plasma  has  been  documented  with  end-on  screak  and  framing 
photographs,  and  radial  Thomson  scattering  scans.  The  formation  of  field 
reversing  layers  (defined  as  (AB^j  > l.OB^)  with  a persistence  time  in 
excess  of  3 been  evidenced  by  local  magnetic  probes  and  several 

diamagnetic  loops.  Typical  layer  cross-section  length-to-width  ratios 
are  equal  to  Aq  cm/2.0  cm  * 20.  That  this  field  configuration  is 
induced  by  plasma  currents,  rather  than  trapped  beam  electrons,  is 
Indicated  by  the  absence  of  hard  x-ray  emission  after  beam  passage. 


The  external  magnetic  field  measurenents  and  diamagnetic  loop  data  are 
fitted  to  an  equilibrium  model  to  estimate  the  radial  profile  of  the 
fields  and  plasma  pressure.  End-on  framing  photograph  neasuremants  of 


the  plasma  radial  dimensions,  Thomson  scattering  measurements  of  nkT 

on  axis  and  2.^  cm  off  axis,  and  CO2  laser  interferometer  determinations 

of  'n.dl,  are  found  to  be  in  agreement  with  this  derived  profile. 

After  beam  passage,  remains  relatively  constant  for  5 .^sec  until 

Bg  becomes  small,  at  which  point  the  annular  plasma  collapses  radially. 

The  measured  1/e  decay  time  of  U.  ^sec,  and  concurrent  two-fold  increase 

in  electron  temperature  (as  determined  by  H line-to-continuum  measure- 

d 

meat s> and  verified  at  one  instant  by  Thomson  scattering)  are  consistent 
with  classical  processes . However,  the  axial  confinement  properties  of 
this  system  cannot  be  inferred  from  the  present  experiment  as  the  L/R 
time  of  the  configuration  is  comparable  to  the  time  for  plasma  particles 
to  be  lost  out  the  ends  of  the  system.  Since  the  ratio  of  L/R  time  to 
plasma  end  loss  time  increases  at  least  as  fast  as  T^  (assuming 
classical  processes),  it  is  necessary  to  increase  the  electron 
temperature  (for  example,  by  an  increase  in  beam  current)  to  conclusively 
assess  the  long-term  confinement  properties  of  this  system. 

II.  Experimental  Apparatus 

The  experimental  apparatus  is  shown  in  Fig.  3.  A rotating 
relativistic  electron  beam  is  produced  by  injecting  an  annular  beam  from 
the  TRITON  accelerator^  (1.2  MeV,  6O  kA,  cm  diameter,  60  nsec  duration) 
through  a non-adiabatic  magnetic  iron  cusp with  a magnetic  transition 


half-width  of  12  eta.  The  applied  axial  magnetic  field,  B^,  approaches 
uniformity  at  distances  greater  than  10  cm  from  each  side  of  the  cusp, 
and  has  a magnitude  of  1.3  kG.  The  stainless-steel  drift  tube  is  1^.6 
cm  in  diameter  and  terminated  30  cm  from  the  cusp  with  a Faraday  cup/ 
calorimeter  assembly.  At  the  usual  fill  pressure  of  33  mTorr  H2  (neutral 
atom  density  of  3*3  X 10^^  cm*^) , net  current  generation  was  found  to  be 
maximized.  Measurements  with  a miniature  rotatable  Faraday  cup  (similar 
in  design  to  that  described  by  Greenspan*^)  and  observations  with  thin 
polyethylene  damage  rods  (similar  to  those  employed  by  Roberson^^)  show 
the  beam  has  approximately  equal  azimuthal  and  axial  velocity  components 
and  has  a typical  radius  of  3.3  cm.  These  diagnostics,  as  well  as  damage 
on  the  main  Faraday  cup  foil,  also  showed  the  axis  of  the  beam-plasma 
system  suffered  shot-to-shot  radial  variations  of  as  much  as  ± 2 cm. 

These  variations  are  thought  to  be  due  to  motion  of  the  beam  in  the 
15  cm  long,  35  mTorr  drift  region  between  diode  and  cusp. 

III.  Experimental  Results 

End-on  streak  and  framing  photographs  were  taken  through  apertures 
mounted  at  the  end  of  the  drift  tube  opposite  the  diode.  (In  these 
experiments,  the  conducting  end  wall  was  defined  by  a transparent  stain- 
less steel  grid,  rather  than  the  Faraday  cup/calorimeter  assembly.)  As 
can  be  seen  from  Fig.  U,  the  plasma  has  an  annular^prof ile  with  radial 
dimensions  comparable  to  those  of  the  rotating  beam  (outer  radius  3.5 
cm,  inner  radius  2.0  cm)  and  is  situated  well  away  from  the  drift  tube 
wall  (radius  7.3  cm).  Mote,  the  plasma  does  not  appear  to  undergo  any 
gross  radial  motion  as  it  decays;  rather,  the  annular  region  appears 
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to  fill  in  and  then  fade  away.  Streak  photographs  of  extended  duration 
show  the  plasma  light  lasts  as  long  as  100  asec,  whereas  the  annular 
profile  is  maintained  typically  for  k to  6 usec.  Supplementary  side- 
on  streak  photographs  (not  shown)  taken  through  a viewing  port  located 
30  cm  from  the  cusp  also  demonstrate  the  plasma  to  be  annular  and  exhibit 
no  gross  radial  motion. 

Thomson  scattering^*^  measurements  at  2.5  cm  off  axis  and  UOO  nsec 
after  beam  injection  show  the  plasma  electrons  to  have,  typically, 
densities  of  5 X 10~^  cm'^  and  a Maxwellian  velocity  distribution 
corresponding  to  temperature  of  5 cV.  At  these  low  temperatures,  the 
plasma  is  expected  to  be  located  only  In  the  light-emitting  regions,  and 
radial  scans  with  Thomson  scattering  is  consistent  with  the  end-on 
photographs  in  Fig.  h.  Because  of  the  shot-to>shoc  variations  of  the 
electron  beam  and  resulting  plasma,  detailed  measurements  of  the  radial 
profile  could  not  be  made.  However,  a series  of  kO  identical  shots,  20 
with  the  scattering  system  probing  on  the  drift  tube  axis  (r  " 0),  and 
20  at  2.5  cm  off  axis  show  (errors  indicate  RMS  scatter); 

-rr:  (r  * 0)  - 7.2  ± 5.2  x 10^®  eV/cm^,  (3a) 

nKl 

ex 


and 


2.5  cm) 


17.1  ± 15.0  X 10^®  eV/cm^. 


(3b) 


(The  higher  nkT 

ei 


at  2.5  cm  off  axis  was  manifested  by  an  increase  in 


boch  eleccroa  temperature  and  density.)  3ata  were  taken  between  '*00  and 


300  nsec  after  beam  injection.  Note  that  the  observation  that  the 

average  nkT^  on  axis  is  greater  than  zero  does  not  imply  nkl^  is 

greater  than  zero  at  the  axis  of  the  rotating  beam:  the  shot-to-shot 

radial  variations  of  the  beam-plasma  system  and  the  finite  radial  extent 

(k  mm)  of  the  scattering  volume  sometimes  bri  g the  confinement  region 

into  the  scattering  volume,  which  is  fixed  in  the  drift  tube  frame. 

(In  fact,  10^  of  the  shots  on  axis  show  nkT  is  less  than 

ej. 

5 X lO"^"*  eV/cm^,  that  is  essentially  zero.) 

Field  reversal  is  inferred  from  the  output  of  the  diamagnetic  loop 
(11.0  cm  in  diameter),  which  measures  the  change  in  axial  magnetic  flux, 


Ac?  (Fig.  5a).  As  the  measured  plasma  pressure,  nkT  , can  account  for 

only  5^  of  the  observed  Acp  (assuming  nkT^  as  in  Eq.  3b,  and  the 

annular  radial  profile  depicted  in  Fig.  4b),  Che  bulk  of  the  diamagnetic 

signal  is  due  to  the  induced  axial  magnetic  field,  <^^2*  ^ploying  a 

simple  flux-conserving  model, in  which  Che  aximuthal  plasma  current 

layer  has  zero  thickness,  is  located  at  a radius  r^,  and  is  inside  a 

flux  conserving  region  of  radius  r^,  the  voltage  induced  on  an  N-turn 

loop  of  radius  r,  (such  that  r > r,  > r ) terminated  with  passive 
L w L P 

integrator  network  of  time  constant  RC,  is  (MRS  units): 


Nrr  AB  r^d-rf/r^) 
z w L w 

RC  (r^/r=-l)' 

w p 


Using  Eq.  ^4-,  the  diamagnetic  signal  depicted  in  Fig.  is  sufficient  to 
account  for  full  cancellation  of  the  applied  field  on  axis  (AB^  » -1.0  B^) 
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for  approxlmacely  6 tisec  (assuming  « 3.5  cm  as  indicated  by  the  end- 

on  photographs).  At  peak  diamagnetic  signal,  > -3.3  kG  or 
approximately  -1.3  The  loss  of  field  reversal  is  concurrent  with  the 

observed  collapse  of  the  annular  plasma  in  Fig.  ^b. 

The  persistence  time  of  the  axial  plasma  current,  as  indicated  by 
the  output  from  a B^  probe  located  6.3  cm  off  axis,  is  shown  in  Fig.  5b. 
A characteristic  feature  of  these  experiments  Is  that  no  appreciable 
decay  of  AB^  Is  seen  until  B^  goes  to  zero.  This  result  is  discussed 
in  Section  V.  These  measurements  show  that  I ^ “ ^5  kA, implying  the 

plasma  currents  induced  by  the  entrance  of  the  beam  are  smell  compared 
to  the  beam  current. 

The  observed  Induced  change  in  axial  magnetic  field  is  consistent 
with  other  measurements.  Using  Eq.  1,  with  v /v„  ■ 1 (determined  from 

Z U 

damage  rod  and  miniature  Faraday  cup  studies),  * ^5  and 

r =3.5  cm,  the  Induced  change  is  calculated  to  be  AB  ■ -2.9  kG  or 
p z 

approximately  -1.6  3^. 

To  satisfy  the  skeptics,  existence  of  field  reversal  was  verified 
directly  by  measurements  with  a miniature  (.3  Cui  diameter)  local  axial 
magnetic  probe.  The  probe  was  insulated  in  a quartz  cube,  positioned 
on  drift  Cube  axis  at  the  same  axial  position  as  the  diamagnetic  loop, 
and  cantilevered  from  a rigid  end  plate  located  1 meter  downstream.  It 
is  generally  desirable  to  obtain  measurements  of  AB^  with  such  an  on- 
axls  probe,  as  opposed  to  a diamagnetic  loop,  since  Che  former  yields 
AB^  directly,  whereas  the  latter  only  if  a plasma  radius  is  known.  The 
probe  output  agrees  well  with  the  diamagnetic  loop,  however,  as  seen  in 
Fig.  6,  lower  field  reversals  (AB^  approximately  equal  to  -1.3B^)  and 
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faster  decay  times  (approximately  2 txsec)  were  observed  when  this  probe  was 
installed.  This  effect  is  attributable  to  the  shot-to-shot  radial 
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variations,  which  would  allow  this  probe  to  intercept  a portion  of  the 
rotating  beam  and  plasma  currents.  (There  is  evidence  in  the  reversed 
field  theta  pinch  experiment  of  Linford*^  that  an  axial  probe  does  not 
have  to  physically  intercept  the  rotating  currents  in  order  to  have  a 
perturbing  effect.)  For  this  reason  the  on-axis  probe  was  not  installed 
during  most  of  this  investigation.  A local  axial  magnetic  probe  was  also 
positioned  at  6.3  cm  radius  (1  cm  inside  the  drift  tube  wall)  to  determine 

outside  the  field  reversing  layer.  Such  a probe  did  not  perturb 
the  plasma  and  indicated  * + 300  G or  approximately  + .17  B^, 
consistent  with  flux  conservation  inside  the  drift  tube. 

Three  concentric  diamagnetic  loops  with  identical  axial  positions 
and  RC  networks  were  employed  to  determine  the  flux  excluding  area  of 
the  current  layer.  If  this  area  is  totally  enclosed  by  the  smallest 
of  these  diamagnetic  loops,  then  by  Eq.  the  signal  voltages  from  all 
three  loops  are  proportional  to  (1  - rf/r^) . In  Table  1,  normalized 
outputs  of  the  three  diamagnetic  loops  are  compared  with  the  calculated 
values.  Agreement  between  calculation  and  experiment  for  these  several 
shots  indicates  the  flux  excluding  area  is  always  inside  the  innermost 
loop  (radius  '^.9  cm).  Allowing  for  the  maximum  radial  variation  of 
- 2 cm,  these  results  suggest  the  layer  radius  cannot  exceed  3.9  cm,  in 
keeping  with  the  previously  estimated  value  of  3.5  cm. 

That  these  field  reversing  layers  are  induced  primarily  by  plasma 
currents  rather  than  trapped  relativistic  beam  electrons  can  be  ascertained 
from  hard  x-ray  studies.  The  emission  from  a 1 mm  diameter  tungsten 
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wire  placed  across  the  drift  tube  at  a position  30  cm  from  the  cusp  is 
shovm  in  Fig.  7>  Appropriate  filters  were  used  to  discriminate  against 
x-rays  with  energies  less  than  70  keV.  The  presence  of  the  tungsten 
target  does  not  affect  either  the  layer  lifetime  or  strength.  As  the 
target  absorbs  all  of  the  beam  electrons  within  1 p.sec,  this  result 
suggests  trapped  beam  electrons  play  a minor  role  in  layer  persistence. 
Furthermore,  comparison  with  Figs.  and  5^  shows  the  emission  exhibits 
neither  the  lifetime  nor  the  temporal  dependence  of  either  the  azimuthal 
or  axial  currents. 

CO2  laser  interferograms^ ' , used  to  determine  the  plasma  line  density 
along  a major  diameter  at  an  axial  position  of  30  cm,  are  shown  on  two 
different  time  scales  in  Fig.  0.  Due  to  limitations  of  the  detector 
response,  the  initial  rise  of  the  plasma  on  the  beam  time  scale  cannot 
be  seen.  After  a comparatively  featureless  period  of  10  ^sec,  during 
which  the  onset  of  "partial  fringes"  is  observed,  the  interferometer 
exhibits  a normal  decay,  as  manifested  by  several  full  amplitude  fringes. 

It  was  ascertained  that  the  quiescent  phase  corresponds  to  a period  when 
the  plasma  refracts  the  COa  laser  beam  out  of  the  diagnostic  arm  path. 

(The  mechanism  for  this  refraction  is  not  v/ell  understood.)  In  keeping 
with  the  end-on  photography  (Fig.  u)  the  partial  fringes  may  correspond 
to  radial  contraction,  but  not  gross  radial  motion  (i.e.,  variation  in 
radial  position)  of  the  plasma.  Assuming,  then,  chat  only  Che  full 
amplitude  fringes  correspond  to  a decay  in  density,  the  resulting  plasma 
line  density,  averaged  over  20  shots,  is  inferred  to  be; 

’n.di  » 66  i 13  X 10^^  electrons/cm^ . 

Note  from  Fig.  3b  plasma  is  present  for  some  100  usec  after  beam  injection. 


a result  consistent  with  extended  end-on  streak  photographs. 

IV.  Comparison  with  Equilibrium  Model 

A toroidal  equilibrium  model  has  been  employed  to  determine  the 

equilibrium  configuration  consistent  with  the  experimental  data.  Input 

to  the  model  include  the  flux  change  measured  by  the  diamagnetic  loop, 

and  the  measured  AB  and  B.  fields  determined  by  local  magnetic 

z 9 

probes  placed  just  inside  the  drift  tube  wall.  The  model  output  is 
found  to  be  in  agreement  with  the  experimentally  determined  plasma 
radial  energy  density  profile. 

In  this  model,  the  plasma  pressure  is  assumed  to  be  balanced  by  the 
Lorentz  force  at  all  points: 


j X B = (5) 

A cylindrical  (r,9,z)  coordinate  system  with  e.  the  unit  vector 

9 

in  the  toroidal  (azimuthal)  direction  (see  Fig.  9)  is  chosen.  The 

plasma  boundary,  where  the  plasma  pressure  is  zero,  is  assumed  to  be 

a rectangle,  whose  center  is  located  at  r = r and  z = 0,  The  cross- 

P 

section  dimensions  are  such  that  -b  s:  z s b and  r - asr<r  + a. 

P P 

Consistent  with  the  assumption  of  axisymmetry  (S/d9  • 0) , the 
magnetic  field  is  expressed  in  terms  of  the  poloidal  (axial)  flux 
function  \|((r,z): 


mi 

r 


0 r 


(6) 
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where  1(7)  Is  a function  of  the  poloidal  flux.  With  the  field 
expressed  in  this  form  and  the  use  of  Maxwell's  equations,  force 
balance  (Eq.  5)  can  be  written  as  a partial  differential  equation  for 
■Jt  (r,z) 


_ ^ i 

3r  r dr  ^^2 


i^TTrS  ^ - I . 


(7) 


To  proceed  further,  it  is  necessary  to  know  P ■ P(7)  and 
I = I(i(r).  Following  Maschke*^,  a quadratic  dependence  for  P and  I 
is  assumed  with  normalization  such  that 


P - Pi  t^  ^ 

o 


(8) 


where  Pi,  and  M are  constants.  In  addition,  the  variation  in 

the  z direction  is  assumed  to  be  harmonic,  i.e.: 


,|f  (r,z)  = t (r)  cos  (kz)  . 


(9) 


The  boundary  condition  that  P=0  at  z=ib  implies  that 
kb  = (2m  + l)(Tr/2).  In  this  calculation  m is  taken  to  be  zero  and,  to 
be  consistent  with  axial  diamagnetic  loop  data  to  be  presented  in  Section 
V,  b =•  X/2  * 20  cm.  Eq.  (7)  can  be  put  in  a more  convenient  form  by 
changing  the  radial  variable  to  p = j2v:V-^r^  . Using  Eqs.  (3)  and  (9), 
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it  is  found  that: 


Here  is  a characteristic  eigenvalue,  where: 


(10) 


k^-M 


(11) 


Although  this  equation  has  analytic  solutions^®,  it  is  more  convenient 

to  solve  it  by  numerical  integration.  Integration  is  carried  out  using 

the  experimentally  determined  values  of  AB  and  B.  at  6.3  cm  off 

z 0 

axis  and  the  change  in  flux  inside  r = 5*5  cm  (i.e.,  the  diamagnetic 
loop  signal)  to  infer  the  radial  position  (at  z = 0)  of  the  outermost 
closed  flux  surface,  then  integrated 

radially  inward  to  r = - a. 

In  Figure  10,  the  calculated  curves  of  the  radial  variation  of 

nkT  (defined  as  0.5  P)  and  B are  shown  as  the  solid  and  broken 
ex  z 

curves,  respectively.  Best  fit  to  the  observed  magnetic  field  data  is 
found  for  = 2.5  cm,  a = 1.75  cm  and  T)  = 0. 

As  sumnarized  in  Table  2,  the  predictions  of  this  model  are  found 
to  be  in  agreement  with  the  experimental  observations  at  t “ ^*-00  to 
300  nsec.  (The  plasma  characteristics  were  constant  over  this  time.) 
The  expected  average  values  for  ^ fixed  radial  position 

are  estimated  from  the  predicted  electron  energy  density  radial  profile 
(solid  curve.  Fig.  10),  and  a C&P  random  number  generator  (to  simulate 


Che  shoc-co-shoc  radial  variaClons).  Based  on  damage  pattern  and  streak 
photographs,  it  is  assumed  that  the  plasma  is  formed  away  from  Che  axis 
with  displacement  6r,  Che  average  of  which  is  (5r)  > 1.0  cm  and  Che 
maximum  of  which  is  l6rl  » 2.0  cm.  Using  20  random  numbers  for  5r, 
normalized  such  that  (6r)  = 1,  the  expected  pressure  at  r - (r)  + 6r  is 
calculated  for  (r)  * 0 and  (r)  ■ 2.5  cm.  These  results  are  in  agree- 
ment with  the  measured  values,  as  are  the  predictions  of  Che  minimum  and 

maximum  energy  densities.  Furthermore,  the  line  density,  I'n.dl, 

u 

(derived  assuming  n^  and  have  Che  same  radial  dependence)  is 

consistent  with  COs  laser  interferometry  data. 

The  magnetic  field  profile  (broken  curve  in  Fig.  10)  predicts 

AB  “ -3.9  kG,  or  -2.17  B , a somewhat  stronger  reversal  than  that 
z o 

calculated  from  the  simple  flux  conservation  model  in  Eq.  4.  The 
discrepancy  is  due  to  Che  assumption  of  an  infintesimally  Chin  current 
layer  in  Che  flux  conservation  model  and  Che  value  from  Che  equilibrium 
model  is  more  reliable. 

The  reversed  field  configuration  does  not  exhibit  any  gross  unstable 
motions.  However,  after  formation,  the  configuration  undergoes  small 
amplitude  (-  10^)  radial  oscillations,  as  manifested  on  Che  diamagnetic 
loop  output.  The  oscillation  period  is  consistent  with  that  predicted 
by  the  theory  of  Chu,  et  al.^°,  in  which  it  was  shown  that  a rotating 
beam  can  excite  the  magneCosonic  mode  of  the  plasma.  The  period  of 
this  mode  is  predicted  to  be: 

T = (2a/V  ) fl  + (nc)2/(2<i,  , (12) 

m L J 
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where  lu  . is  Che  ion  plasma  frequency,  V ■ V <1  +3)®  (V.  Is  Che 
pi  m A A 

Alfven  velocicy,  9 ■ 2uQnkT^/B^) , and  I and  r are  respeccively  che 

lengch  and  radius  of  che  beam.  Wlch  che  parameCers  of  chis  experimenc, 

3 » .5,  » 5 X 10‘®  cm*^,  a • 3.5  cm,  2*4-0  cm,  and  ■ 1.2  kG  (ac 

che  plasma/vacuum  incerface,  see  Fig.  10),  ic  is  found  cha^  t ■ 1.6  usee, 

n 

in  good  agreemenc  wich  Che  observed  period  of  approxlmacely  2 usee. 
Calculacions  IndicaCe  che  radial  excenc  of  chese  oscilladons  is  on  che 
order  or  - .5  cm.  Of  imporcance  here  is  chac  chese  oscilladons  do  noc 
grow  in  dme,  suggesdng  chaC  che  syscem  has  "found"  a scable  equilibrium. 

The  experimencal  observadons , Chen,  confirm  Chac  an  equilibrium 
elongaCed  cross-secCion  Coroidal  configuraCion  has  been  produced.  The 
resulcing  coroidal  cross-secdon  is  non-circular  wiCh  a lengch- co-widch 
ratio,  b/a,  approximately  equal  to  20.  The  z = 0 (midplane)  cross- 
section  experimental  average  value  of  3 = < 2 UQrik(T^  + T^^)  /B^  > is 
approximately  .5.  Unlike  a Tokamak-type  configuration,  the  poloidal 
(axial)  field  on  Che  outer  closed  flux  surface  is  larger  chan  the  toroidal 
(azimuthal)  field.  Nevertheless,  che  small  aspect  ratio  (r^/a  on  the 
order  of  I.5)  and  large  length-to-width  ratio  yield  an  effective  value 
for  q of  4 on  the  outer  flux  surface,  where 

q = 4(b  + a)  B,/2rrr^B,  . (13) 

3 p Z 

Thus,  one  would  expect  this  configuration  to  be  magnetohydrodynamically 
scable . 
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V.  The  Decay  of  the  Configuration 

After  the  beam  has  left,  the  induced  magnetic  configuration  decays 

as  a result  of  the  resistive  dissipation  of  the  plasma  currents.  The 

signals  from  the  diamagnetic  loop  (which  measures  external  to  the 

loop)  and  the  probe  (Fig.  3)  not:  decay  exponentially  because  the 

changing  magnetic  fields  cause  a change  in  the  plasma  equilibrium.  The 

decay  of  j decreases  b|,  changing  the  force  balance  to  move  the 

plasma  outward,  whereas  the  decay  of  will  cause  B^  on  axis  to 

decrease  faster  than  B^  outside  the  layer,  and  thus  move  the  plasma 

inward.  Thus  it  is  reasonable  to  expect  a net  outward  motion  of  the 

plasma  which  compresses  the  (decaying)  axial  return  flux  against  the 

wall,  keeping  B outside  the  loop  approximately  constant  until  B.  has 
z 9 

become  small,  at  which  point  the  plasma  column  contracts.  This 
explanation  is  consistent  with  the  framing  photographs  (Fig.  kb),  which 
show  an  increase  in  the  mean  plasma  radius  in  the  first  k usee  and 
subsequent  contraction;  and  a simple  pressure  balance  model,  which 
correctly  predicts  the  temporal  behavior  of  the  diamagnetic  loop  signal 
shown  in  Fig.  ’^a.. 

The  observed  late  time  1/e  decay  times  for  both  B,  and  AB  of 

0 z 

approximately  k usee  are  consistent  with  classical  resistive  decay  of 
the  currents,  assuming  the  electron  temperature  of  5 eV  measured  by 
Thomson  scattering.  However,  at  these  low  temperatures,  the  L/R  time 
of  the  configuration  is  also  comparable  to  t^,  the  time  for 

plasma  to  free  stream  out  of  the  ends  of  the  apparatus  (given  by  x/2c 


where  I is  the  length  and  c^  the  sound  speed) , assuming  no  axial 
confinement.  Hence,  although  it  is  believed  the  field  lines  are  closed 
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at  the  ends,  this  set  of  experiments  does  not  provide  direct  evidence 
of  plasma  containment.  Since 


3/2  -1/: 

Tjr  « T /T 
D s e e 


(li^) 


to  verify  axial  containment  in  the  present  experiment  requires  T to 

e 

exceed  30  eV. 

As  is  to  be  expected,  the  dissipation  of  magnetic  energy  leads  to  an 

increase  in  plasma  temperature,  as  indicated  by  the  data  in  Fig.  11, 

which  shows  the  results  of  H line-to-continuum  measurements.  Light 

3 

was  collected  along  a chord  2.5  cm  off  the  drift  tube  axis,  and  wave- 

0 

length  resolved  with  three  27  A wide  channels,  one  at  the  H line 
0 O 3 

(^<•861  A)  and  the  other  two  at  4a6l  i 27  A.  As  seen  in  Fig.  11,  the 

electron  temperatures  indicated  by  this  diagnostic  are  in  agreement  with 

those  predicted  by  Thomson  scattering.  The  most  salient  feature,  however, 

is  the  30^-  increase  in  T over  the  layer  lifetime.  The  amount  of 

e 

electron  heating  that  can  take  place  through  the  induced  current  decay 

can  be  calculated  assuming  classical  resistivity.  Using  Che  measured 

initial  electron  temperature  of  5 sV  at  full  ionization,  calculations 

of  the  total  energy  deposited  into  the  plasma,  i.e.,  J^TItV,  where 

is  the  axial  current,  T)  the  classical  resistivity,  t the  appropriate 

time,  and  V the  appropriate  volume,  yield  a AE/pair  = T + T.  of 

e 1 

approximately  12  eV.  Assuming  a classical  electron-ion  relaxation 
time  (T^  = T^  on  these  time  scales)  yields  a final  electron  temperature 
of  approximately  11  eV,  consistent  with  that  observed  in  Fig.  11. 

It  is  not  surprising  that  the  resistivity  observed  in  these  experiments 
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is  cLasslcaL,  as  aC  these  low  temperatures,  single  part'cle  Coulomb 
collisions  dominate.  For  the  classical  collision  frequency  to  be  less 
chan  Che  lower  hybrid  frequency,  for  example,  Che  electron  temperature 
must  exceed  55  (The  lower  hybrid  mode^^  is  the  most  applicable 

instability  for  this  geometry  as  it  yields  an  anomalous  cross-field 
resistivity  and  can  readily  be  excited  whenever  the  drift  velocity  exceeds 
the  ion  thermal  velocity.) 

The  layer  length  appears  to  decrease  as  the  system  decays.  The  out- 
puts of  three  diamagnetic  looos  placed  at  various  axial  positions  are 
shown  as  a function  of  time  in  Fig.  12.  At  500  nsec  after  beam  injection, 
the  layer  strength  kO  cm  from  the  cusp  is  approximately  30'^  of  chat  at 
10  cm  from  Che  cusp,  suggesting  a characteristic  length  of  cm.  The 
signals  from  Che  loops  at  the  ends  of  the  layer  decrease  faster  than 
chose  from  Che  one  near  the  center.  This  observation  can  be  explained 
by  modeling  Che  azimuthal  currents  as  a coil  which  tends  Co  maximize  its 
self  inductance,  i.e.  contract  axially.  This  tendency  is  balanced  by 
Che  pressure  of  Che  confined  plasma,  which  maintains  Che  elongated 
geometry.  (This  may  explain  why  the  field  reversing  layers  in  this 
experiment  do  not  coalesce  into  a "bicycle  tire"  geometry  characteristic 
of  field  reversing  relativistic  electron  rings^^,  which  do  not  contain 
a warm  plasma.)  Therefore,  as  plasma  pressure  is  lost,  Che  layer  length 
should  decrease,  as  observed. 

VI.  Energy  Accounting 

An  attempt  was  made  to  ascertain  the  partitioning  of  the  rotating 
beam  energy  after  passage  through  its  self-made  plasma.  The  energy 


accounting  is  not  yet  complete,  as  can  be  seen  in  Table  3.  These  numbers 
were  calculated  at  ^00  nsec,  when  the  system  is  thought  to  be  in  quasi- 
steady  state.  A layer  o£  inner  and  outer  radii  2.0  and  U..0  cm, 
respectively,  and  axial  length  ^0  cm  was  used  in  these  calculations. 

Of  the  original  1^20J  in  the  beam,  40^  is  lost  in  passing  through  the 
cusp.  This  low  transfer  efficiency  is  thought  to  be  due  to  the  broad 
cusp  transition  width  and  long  drift  region  between  diode  and  cusp.  Of 
the  injected,  20^  is  deposited  in  the  plasma  and  induced  fields, 

passes  straight  through  the  system  and  38*^  is  unaccounted  for.  One 
of  the  most  important  tasks  of  future  experiments  is  to  account  for  the 
"lost"  energy  and  to  improve  the  fraction  of  beam  energy  that  is  deposited 
in  the  plasma.  However,  using  the  energy  accounting  of  Table  3,  some 
information  may  be  inferred  about  the  rotating  beam-plasma  energy  transfer 
during  the  beam  pulse: 

1.  With  a cathode  radius  of  2 cm  and  an  observed  damage  radius  on 
polyethylene  rods  of  3.5  cm,  radial  expansion  producing  large  amplitude 
Alfven  wave  heating"^  is  a possible  energy  transfer  mechanism.  This  in 
turn  implies  the  ions  can  have  a substantial  temperature.  This  aspect 
will  be  investigated  further  in  a future  experiment. 

2.  High  electron  temperatures  are  observed  on  a time  scale  less 
than  four  times  the  beam  pulse  duration,  indicating  a rapid  beam-to- 
plasma  energy  transfer.  Classical  processes,  such  as  Coulomb  scattering 
of  beam  primary  electrons  or  classical  resistive  dissipation  of  induced 
plasma  return  currents,  cannot  account  for  such  a rapid  energy  transfer. 

In  the  latter  case,  computations  of  ionization  and  heating  by  return 
current  dissipation  have  been  performed  using  a code  described  else- 
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where. It  is  found  that  such  rapid  LonizaCion  and  heating  in  Che 
presence  of  classical  resistivity  alone  would  require  plasma  current 
densities  in  excess  of  I5  kA/cm^,  an  order  of  magnitude  larger  chan 
can  reasonably  be  expected  in  this  experiment.  On  the  ocher  hand, 
turbulent  resistivities,  obtained  from  models  of  Che  e-e,  ion  acoustic, 
or  lower  hybrid  drift  modes,  are  capable,  either  individually  or  in 
combination,  of  producing  the  observed  ionization  levels  and  electron 
temperatures  with  current  densities  of  only  2 kA/cm.*  This  current 
density  is  more  in  keeping  with  chose  inferred  from  Che  present 
experiment  (0.8  - 1.6  kA/cm^) . 

3.  It  should  be  noted  again  chat  Che  beam  deposits  most  of  its 
energy  into  the  induced  magnetic  field  (103  J of  magnetic  energy  vs 
24  J of  plasma  energy)  and  that  dissipation  of  the  induced  currents 
is  a significant  heating  mechanism  in  this  experiment. 
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Table  I 


Comparison  of  calculated  and  observed  diamagnetic  signal  for  the  radially 
resolving  diamagnetic  loops.  Observed  signals  are  taken  from  four  shots, 
at  varying  times  after  electron  beam  injection.  All  data  are  normalized 


Quantity  j 

! 

Method  of 

Experimental 

Theory  ! 

— 

Determination 

Observations 

• 

— — ' t 

1 

Radial 

End-on  Framing 

r;^  • 2 cm 

r:.  • 1.5  cm 

Otmenslons 

Photographs 

r^  » 3.5  cm 

cz  • k cm 

Position  of 

Damage  Rods , 

r • 5 . 5 cm 

r • 3.25  cm 

'^9  Maximum 

Concentric  Loops 

P 

P 

i nkT  (r  - 0) 

1 

Thomson  Scattering 

T.2i5.5Xl0-^eV/cm^ 

d 

6.5il.Txl0^^eV/cm^ 

Table  II 

Sunriary  of  comparison  between  experimentally  determined  and  theoretically 
predicted  quantities,  based  on  an  equilibrium  model  consistent  with  observed 
magnetic  field  data.  Data  taken  between  400  to  600  nsec  after  beam  injection 


nKT  ( r«2 . 5 cm  ) 
ei 


nkT  . (r-0) 
eimin 


nkT  (r“2.5  cml 
eimax  1 


I lT.ltl5.0xl0^^eV/cm^  - 25tl.i4Xl0’‘^eV/cm^ 


b 


<5xl0-‘*eV/cm^ 


' oOxlO^^eV/cm^ 


55^iO*^eV/ctn 


iCO^Interferomecer  06^13/ I0*^^elec/cm^  l60x  lO^^elec/cm 


|dB  on  axis  jDiamagneClc  Loop  -1.6  kC 

I 2 i 


a Average  of  20  shots 
b Average  of  three  lowest 
c Average  of  three  highest 
d 20  "random"  shots 


REVERSED  MAGNETIC  FIELD  CONFIGURATION 
RESULTING  FROM  INDUCED  AXIAL  MAGNETIC  FIELD 


Fig.  1 — Rotating  relativistic  electron  beam  induced  axial 
*jz-net^  azimuthal  plasma  currents  and  the  re- 

sultant azimuthal  (toroidal)  field,  B^,  and  axial  (poloidal) 
field  respectively 


Fig.  2 — Creation  of  reversed- field  configuration  by  a rotat- 
ing relativistic  electron  beam:  (a)  beam  enters  initially  neu- 
tral gas  producing  self-magnetic  fields;  (bl  plasma  is  formed 
and  heated;  (c)  beam  leaves  the  system,  inducing  plasma 
currents  which  persist  in  the  high  conductivity  plasma,  result- 
ing in  further  plasma  heating  as  they  decay.  (For  simplicity, 
the  induced  and  resultant  outside  the  confined  plasma 
have  not  been  shown.) 
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Fig.  4 — End-on  streak  (a)  and  framing 
(b)  photos  of  plasma  column 
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Fig.  5 — Typical  oscilloscope  traces  showing  magnitude  and 
temporal  behavior  of  the  induced  (a)  axial  and  (b)  azimuthal 
fields.  Trace  (a)  is  the  diamagnetic  loop  signal.  The  lower 
dashed  line  indicates  the  diamagnetism  accounted  for  by  the 
electron  energy  density  as  determined  from  Thomson  scatter- 
ing. The  remeiinder  of  the  diamagnetic  signal  is  due  to  rotat- 
ing beam-induced  field  reversing  plasma  current  layers.  The 
right  hand  vertical  scale  gives  the  change  in  axial  field  on  axis, 
assuming  a layer  radius  of  3.5  cm,  in  terms  of  the  applied 
field,  Bjj.  Trace  (b)  is  calibrated  both  in  units  of  the  azi- 
muthal field  (kG)  and  the  axial  current  (kA). 
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Fig-  8 — CO2  laser  interferometer  traces  on  1 /isec/div  (left) 
and  10  jusec/div  (right) 
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Fig.  9 — Geometry  used  in  equilibrium  model 
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Fig.  10  — Electron  energy  density  profile  (solid  curve)  and 
axial  magnetic  field  profile  (broken  curve)  predicted  by  an 
equilibrium  model  fitted  to  the  magnetic  field  measurements. 


36 


nkTe,  x 10  ^ (ev/cm' 


